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Abstract 


The Fractal Flame algorithm is a member of the Iterated Function 
System (IFS) class of fractal algorithms. A two-dimensional IFS creates 
images by plotting the output of a chaotic attractor directly on the image 
plane. The fractal flame algorithm is distinguished by three innovations 
over text-book IFS: non-linear functions, log-density display, and struc- 
tural coloring. In combination with standard techniques of anti-aliasing 
and motion blur the result is striking image variety and quality. 

The guiding principle of the design of the algorithm is to expose and 
preserve as much of the information content of the attractor as possible. 
We found that preserving information maximizes aesthetics. 


1 Overview 


Some examples appear in Figure 1. The paper begins by defining classic, linear 
iterated function systems and hence grounding our notation and terminology. 
The classic formulation is then extended with non-linear variations in Section 3, 
and then further extended with post transforms and final transforms in Sec- 
tion 3.1. Section 4 describes how log-density display works and its importance, 
and Section 5 covers the coloring algorithm. These three sections cover the core 
innovations. Sections 6 and 8 explain other important properties of the imple- 
mentation, and Section 7 shows how to create symmetric flames. The appendix 
is a catalog of the variations including formulas and examples. 


2 Classic Iterated Function Systems 


A two-dimensional Iterated Function System (IFS) is a finite collection of n 
functions F; from В? to R2. The solution of the system is the set S in R2 (and 
hence an image) that is the fixed point of Hutchinson’s recursive set equation 
[3]: 


S= U Е(5) 
1=0 
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Figure 1: Example fractal flame images. The names of these flames are: a) 206, 
b) 191, c) 4000, and d) 29140. These images were selected for their aesthetic 
properties. 
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Figure 2: Sierpinski’s Gasket, a simple IFS. X and у increase towards the lower 
right. The recursive structure is visible as the whole image is made up of three 
copies of itself, one for each function. 


As implemented and popularized by Barnsley [1] the functions F; are linear 
(technically they are affine as each is a two by three matrix capable of expressing 
scale, rotation, translation, and shear): 


Fi(z,y) = (aix + biy + ci, dia + eiy + fi) 


For example, if the functions are 


Fo(z, y) = (5 5) Кі.(а,у) - (==, 2) Р(х,у) = ev ун) 


then the fixed-point 5 is Sierpinski’s Gasket, as seen in Figure 2. 

In order to facilitate the proofs and guarantee convergence of the algorithms, 
the functions are normally constrained to be contractive, that is, to bring points 
closer together. In fact, the normal algorithm works under the much weaker 
condition that the whole system is contractive on average. Useful guarantees 
of this become difficult to provide when the functions are non-linear. Instead 
we recommend using a numerically robust implementation, and simply accept 
that some parameter sets result in better images than others, and some result 
in degenerate images. 

The normal algorithm for solving Юг S is called the chaos game. In pseu- 
docode it is: 


(x, y)= а random point in the bi-unit square 
iterate { 
i = a random integer from 0 to n — 1 inclusive 
(т, y) = Fi(z, y) 
plot(z, у) except during the first 20 iterations 


The bi-unit square are those points where x and y are both in [-1,1]. The 
chaos game works because if (x,y) € 5 then Е.(т,у) € 5 too. Though we start 


out with a random point, because the functions are on average contractive the 
distance between the solution set and the point decreases exponentially. After 
20 iterations with a typical contraction factor of 0.5 the point will be within 
10-6 of the solution, much less than a pixel's width. Every point of the solution 
set will be generated eventually because an infinite string of random symbols 
(the choices for i) contains every finite substring of symbols. This is explained 
in more formally in Section 4.8 of [1]. 

No sufficient number of iterations is given by the algorithm. Because the 
chaos game operates by stochastic sampling, the more iterations one makes the 
closer the result will be to the exact solution. The judgement of how close is 
close enough remains for the user. 

In fractal flames, the number of samples is specified with the more abstract 
parameter quality, or samples per output pixel. That way the image quality 
(in the sense of lack of noise) remains constant in face of changes to the image 
resolution and the camera. 

It is useful to be able to weight the functions so they are not chosen with 
equal frequency in line 3 of the chaos game. We assign a weight, or relative 
probability ш; to each function F;. This allows interpolation between function 
systems with different numbers of functions: the additional functions can be 
phased in by starting their weights at zero. Differently weighted functions are 
also necessary to draw symmetric flames as shown in Section 7. 

Some implementations of the chaos game make each function’s weight pro- 
portional to its contraction factor. When drawing one-bit per pixel images 
this has the advantage of converging faster than equal weighting because it 
avoids redrawing pixels. But in our context with mutiple bits per pixel and 
non-linear functions (where the contraction factor is not constant) this device 
is best avoided. 


3 Variations 


We now generalize this algorithm. The first step is to use a larger class of 
functions than just affine functions. Start by composing a non-linear function 
V; from R? to R? with the affine functions: 


Fi(z,y) = У, (aie + biy + ci, dix + ei + fi) 


We call each such function V; a variation, and each one changes the shape 
and character of the solution in a recognizable way. The initial variations were 
simple remappings of the plane. They were followed by dependent variations, 
where the coefficients of the affine transform define the behaviour of the vari- 
ation. More recently, parametric variations have been introduced, which are 
variations controlled by additional parameters independent of the affine trans- 
form. 

Appendix A documents 49 variations. Variation 0 is the identity function. 
It and six more examples are: 


Vo(z,y) = (mv) linear 
Vi(z,y) = (sinz,siny) sinusoidal 
Vo(z,y) = = “(а,у) spherical 
V3(a,y) = (xsin(r?) — y cos(r?), x cos(r?) + ysin(r?)) swirl 
У(љу) = =: ((æ-—y)(x +y), 22у) horseshoe 
Viz(z,y) = (x+ csin(tan Зу), у + f sin(tan 32)) popcorn 
Voi(m,y) = (sin(piy) - cos(p22), sin(p3x) — cos(pay)) рај 


where 
= a 

An example of a dependent variation is Popcorn, Ут, which is dependent 
on the c and f coefficients of the affine transform. 

The PDJ variation, V24, is an example of a parametric variation. PDJ relies 
on four external parameters (pi, рэ, pa, ра) to fully characterize its behaviour. 

Variations can be further generalized by replacing the integer parameter j 
with a blending vector v;; with one coefficient per variation. Then 


Е; (т, y) ЭЛС Ї biy | Ci, d;z | сіў | fi) 
J 


With this generalization, if are n functions, then there are 87n parameters 
that specify it. The parameters consist of 49n variational coefficients vij, 31n 
parametric coefficients, 6n matrix coefficients a; through f;, and п weights их. 
Significantly fewer parameters are required, however, if unused parameters are 
assumed to be 0. 


3.1 Post Transforms 


To this point, applying a transform to a set of coordinates involves first applying 
an affine transformation to the coordinates, and then applying a non-linear vari- 
ation function to the result of the affine transformation. We further generalize 
this with the addition of a secondary affine transformation, a post transform, to 
be applied after the non-linear function. This provides the ability to change the 
coordinate systems of the variations. If the post transform is 


P,(x,y) = (ona + @ у + Yi, дах + ciy + G) 
then we redefine the F; as follows: 


Е(х,у) = P vig Vj (aie + biy + ci, dix + exy + fi)) 
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3.2 Final Transforms 


We can now introduce the concept of a final transform, which is an additional 
function Рета (2, y) that is always applied regardless of the value of i in the 
iteration loop. The final transform is like a non-linear camera. The result of 
the application of the final transform function is not ‘in the loop’, and there 
can be only one final transform per flame. The final transform can have a post 
transform associated with it. In pseudocode, we now have: 


(x, y)= a random point in the bi-unit square 
iterate { 
i = a random integer from 0 to n — 1 inclusive 
(2, y) = ЕЁ, (т, у) 
(жу, уу) = Ета (2, У) 
plot (27, уг) except during the first 20 iterations 


Adding the post transforms to the earlier parameter count, we now have 93n 
parameters necessary to fully specify n functions. If a final transform is desired, 
then the total number of parameters necessary is 93(n + 1). 


4 Log-Density Display 


The chaos game produces a series of (2, y) points which are plotted on the image 
plane. The collection of these points approximates the solution S to the iterated 
function system. 5 is a subset of the plane, and hence membership is a binary 
function, and the image is therefore black and white, lacking even shades of 
gray. See Figure 3a for an example. 

Information is lost every time we plot a point that has already been plotted. 
A more interesting image can be produced if we render a histogram of the 
chaotic process, that is, increment a counter at each pixel instead of merely 
plotting. These counters can be visualized by mapping them to shades of gray 
or by using a color-map that represents different densities (larger count values 
are more dense) with different colors. A linear mapping of counters to gray 
values results in Figure 3b. 

The result is unsatisfying because of the large dynamic range of densities. 
Like many natural systems, the densities are often distributed according to a 
power law (frequency is proportional to an exponent of the value). Figure 4 has 
two histograms demonstrating this. The densest points are much denser than 
the average density, hence with a linear map most of the image is very dark, 
and information is lost. 

The flame algorithm addresses the problem by using a logarithmic map from 
density to brightness. See Figure 3c for the result. The logarithm allows one 
to visually distinguish between, for example, densities of 3,000 and 5,000 in one 
part of the image and 30 and 50 in another part. 
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Figure 3: Successive refinements of the rendering technique starting with a) 
binary membership, then b) linear, c) logarithmic, d) with color, e) with gamma 
factor, and finally f) with vibrant colors. The parameters to create this image 
are given in Appendix B. 
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Figure 4: Plots showing that the distribution of densities in an IFS follows the 
power law. The density (on the horizontal axis) is the number of hits by the 
system in a pixel, the frequency (on the vertical axis) is the number of pixels 
with that density (or up to the next power of two). The line is from the image 
in Figure 3, the other is the composite of 19 old, favorite systems including 
Figure 3. In each case, after a plateau the graph is nearly a straight line in 
log-log space. This is the definitive characteristic of a power law distribution. 
Each image was computed with 9.2e6 samples on a 900x900 grid. 


The display of high dynamic range images like these by tone mapping is 
studied in computer graphics [4]. The logarithm used here (combined with the 
gamma factor, described below) is just an ad-hoc tone-map. 

This log-density mapping is the source of a 3D illusion. On sight people 
often guess that fractal flames are rendered in 3D, but as just described the 
algorithm works strictly with a 2D buffer. However, where one branch of the 
fractal crosses another, one may appear to occlude the other if their densities 
are different enough because the lesser density is inconsequential in sum. For 
example branches of densities 1000 and 100 might have brightnesses of 30 and 
20. Where they cross the density is 1100, whose brightness is 30.4, which is 
hardly distinguishable from 30. 


5 Coloring 


There is more information to be wrung from the attractor. In particular, which 
parts of the attractor come from which functions? The flame algorithm uses 
color to convey this. The result is a substantial aesthetic improvement. Color 
could be assigned according to the density map and although the result is in- 
creased visibility of the densities relative to grayscale, the internal structure of 
the fractal remains opaque. Furthermore, for animation the eye prefers that the 
color of each part of the attractor remain unchanged over time, otherwise the 
illusion of an object in motion is compromised. The fractal flame algorithm uses 
an original means to accomplish this: adding a third coordinate to the iteration. 

Naturally we want to use a palette or color-map which we define as a function 
from [0,1] to (r,g,b) where г, 5, and b are in [0,1]. А palette is classically specified 
with an array of 256 triples of bytes. 

To achieve this we assign a color c; to each function F; (and a corresponding 
Cfinal if a final transform is present) and add an independent coordinate to the 
chaos game: 


(x, y)= а random point in the bi-unit square 
с а random point in [0,1] 
iterate { 
i = a random integer from 0 to n — 1 inclusive 
(x,y) = Fi(z, у) 
c= (с + с;)/2 
(ту, yf) = Ета (2, у) 
су = (с + Сріпа)/2 
plot (27, уў, ср) except during the first 20 iterations 


This has the important property that the most recently applied function 
makes the largest difference in the color index and also in spatial location. 
Indices make less difference as they recede in time. Hence colors are continuous 
in the final image. 


Color plotting is naturally implemented by keeping three counters per pixel 
instead of one, and adding the current color to the three of them instead of incre- 
menting a single density counter. That is not enough information for proper log- 
density display, however. Taking the logarithm of each channel independently 
grossly alters the colors. Instead one must add a fourth channel of so-called 
alpha (а), or transparency values. 

So then to plot a point the color is added to the three color channels, and 
1 is added to the alpha channel. After all the samples have been accumulated 
each channel is scaled by loga/a. See Figure 3d for the result. The resulting 
alpha values can be output with the image if the file format supports them, or 
they can be used for compositing the fractal with a background immediately, or 
they can be discarded. 


6 The Gamma Factor 


Accurate display of any digital image on a Cathode Ray Tube (CRT) requires 
gamma correction to account for the non-linear response of screen phosphors to 
voltage. Without correction the darker parts of an image appear too dark. If 
the brightness б of a pixel is a value between 0 and 1, the corrected brightness 
is simply: beorrected = 0177 where y normally about 2.2. 

But depending on the specific image, gamma values as large as 4 improve vis- 
ibility of fine structure in the attractor. Large gamma values also substantially 
increase the visible noise, and so require longer rendering times to compensate. 
The parameter is therefor left to the user to set according to taste and cir- 
cumstance. See Figure 3e for the result of applying gamma 4 to our running 
example. 

Because the red, green, and blue phosphors respond independently, gamma 
correction is normally applied to each color channel independently. However 
if the gamma value is unnaturally large this has the effect of washing out the 
colors of the image (it becomes pastel or ghostly off-white). This is because 
saturated colors occur due to large difference between channel values. But 
gamma correction boosts any small values towards one, leaving less difference, 
and hence less saturation. 

While one may desire this effect, to preserve bright colors the gamma cor- 
rection may be applied the same way the logarithm is: by calculating a scale 
factor on the alpha channel, and then applying it to the three color channels. 
The name of the parameter that selects this is called vibrancy, and it can take 
any value from 0, meaning to apply gamma to channels independently, to 1, 
meaning to apply gamma from the alpha channel to each channel. 


7 Symmetry 


The human mind responds to symmetric designs at a fundamental level. When 
the matrix coefficients are chosen at random, the chances of a symmetric design 
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appearing are vanishingly small. We can easily inject such functions intention- 
ally, however. The result appears in Figure 5. 

There are two kinds of symmetries: rotational and dihedral. First we cover 
rotations. Adding a function to the system that rotates by 180 degrees makes a 
2-way rotational symmetry appear. The weight of this function should be equal 
to the sum of the weights of all other functions in the system. That way half 
of the jumps in the chaos game are between the two halves, and hence the two 
halves will have equal density. If the rotation function is given the same weight 
as the other functions, then one of the two halves will only be a shadow of the 
other. 

Adding a rotation by 120 degrees does make a 3-way symmetry appear. 
The three branches do not have equal density however, and no weighting would 
balance them. That’s because in order to get into the 240-degree branch in the 
chaos game, one has to pick the rotational function twice in a row, which is only 
25% probable, but the 120-degree branch is 50% probable. Instead one must 
introduce two transforms, one by 120 degrees and one by 240, and give them 
both weight equal to the sum of the others. Then all three branches will have 
the same probability. In general, to produce n-way symmetry, n-1 additional 
transforms are necessary to balance the densities. 

A dihedral symmetry is created by adding a function that inverts the x coor- 
dinate only (producing bilateral symmetry). Again it is given weight equal to the 
sum of all the other weights combined. Combining this function with rotation 
functions gives all the dihedral symmetries. Dihedral symmetries are named 
with negative integers, so the simple bilateral symmetry is -1, and snowflake 
symmetry is -6. This just follows the isomorphism between multiplication on 
integers and composition of symmetries. 

It is also possible to introduce symmetries by modifying the chaos game to 
support them directly instead of adding symmetric functions. For example, one 
can just add an integer symmetry parameter, and then after picking a function 
at random, pick a rotation at random. But then how to interpolate between 
symmetries without discontinuities is problematic. 

Getting good colors with the symmetry functions requires special treatment. 
The problem is the symmetry functions can bring a point back onto itself after 
two or more applications. If the color is modified by these functions and then 
plotted at its original location, different colors get averaged, and the image loses 
color diversity. The solution is to not change the color coordinate when applying 
symmetry functions. This also makes the colors symmetric as well as the shape. 


8 Filtering 


Aliasing in spatial and temporal directions is visually disturbing and also indi- 
cates information loss because one cannot tell if an artifact in the image is an 
original or an alias. With anti-aliasing, there is no ambiguity. 

The chaos game lends itself to anti-aliasing. Consider spatial aliasing first, 
that is the elimination of the jaggie edges. The normal technique is to draw 
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Figure 5: Examples of symmetry. Image d shows how the colors wash out 
without special treatment of the color coordinate for symmetry transforms. 


the desired image at high resolution, and then filter it down to display resolu- 
tion. This is known as supersampling, and its cost is linear in time and memory 
(though the sampling is normally applied in two dimensions, so 3 by 3 super- 
sampling means 9x). With the chaos game however we can achieve this effect 
by just increasing the number of buckets used in the histogram without increas- 
ing the number of iterations. There is a small, sublinear cost in time though 
because the time spent filtering is significant, and the increased memory usage 
also means increased cache misses during iteration. The effect on visual quality, 
however, is dramatic. Gamma correction should be done at this filtering step, 
when the maximum number of bits of precision is still available. 

Despite this filtering, the logarithm and gamma factor may cause low-density 
parts of an image to appear dotted or noisy. A wider filter would solve this, but 
at the expense of making the well-sampled parts of the image blurry. This can 
be addressed with a form of Density Estimation [5]. We have implemented a 
dynamic filter where a blur kernel of width inversely proportional to the density 
of points in the histogram is applied to the samples [6]. The blur kernel is scaled 
based on the supersampling level, and is applied after the log density scaling. 
This variable width filter allows higher density areas to remain in focus, while 
lower density areas are significantly smoothed. Figure 6 illustrates the effect of 
density estimation. 
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Figure 6: Demonstration of density estimation: (a) а low-resolution, low-quality, 
zoomed image rendered without density estimation, and (b) with density esti- 
mation, and (c) a high quality render at high resolution with density estimation. 


9 Motion Blur 


Motion blur, or temporal anti-aliasing, is not so easy to do correctly. Super- 
sampling can be achieved by varying the parameters over time while running 
the chaos game. That is, if 5M samples are used to draw the frame at time t, 
to instead use 1M samples at time t-0.5, 1M samples at t-0.25, 1M at t, 1M at 
1--0.25 and 1M at ++0.5, all the while accumulating into the same buffer. That 
would be 5x supersampling for free! With linear density display, this would 
work exactly. 

The non-linearity of the logarithm complicates things. Consider a pixel 
with density 8. Assume for simplicity the logarithm is base 2, so its assigned 
brightness is 3. Now put the fractal in motion so it blurs across two pixels, each 
should get brightness 1.5. But if the motion takes place before the logarithm 
then the 8 samples would be divided into two pixels of density 4 each, whose 
logarithm is 2, not 1.5. Objects in motion would appear unnaturally bright. 

A proper solution requires the use of an extra buffer: the first buffer is 
linear and accumulates the histogram. After each temporal sample, take the 
logarithm of this buffer and accumulate it into the second one, applying the 
density estimation filter in the process. After all samples are completed, the 
second buffer is filtered down into the final image. 

The drawback of this approach, however, is the computational effort required 
to apply the density estimation filter repeatedly to the linear buffer; multiple 
applications of the filter can easily double the rendering time. 

After experiments with a single buffer yielded acceptable results, we decided 
to default to single buffer renders, while retaining the option of using the extra 
buffer. 


9.1 Directional Motion Blur 


Uniformly distributing the samples among time steps works well for an animated 
series of frames, but the illusion of motion of an individual frame animation may 
be improved by providing a sense of direction. Early attempts at directional 
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Figure 7: (a) motion blur and (b) directional motion blur. 


motion varied the number of samples used at each time step, but since the 
density estimation filter was more aggressive during less dense time steps, the 
result appeared unnatural. 

Instead, the color of points accumulated during earlier time steps are scaled 
in intensity, with the scaling constant approaching 1.0 as the time steps progress. 
A rendering parameter can be supplied to the renderer to change how aggres- 
sively the blending varies throughout the time steps. See the effects of directional 
motion blur in Figure 7. 


10 History and Acknowledgements 


In 1987 at Brown University Bill Poirier showed Scott Draves what he called 
‘Recursive Pictures’ a kind of two-dimensional IFS. Poirier used a formulation 
that included perspective transforms, but lacked a software implementation. In 
response Draves created the first of many IFS algorithms. It was written in 
Postscript and ran on the Laserwriter, producing high resolution line drawings. 
Draves reimplemented this idea in a variety of ways over the three years he 
spent in the Computer Graphics Research Group at while still at Brown. 

The first implementation to include all three definitive characteristics of 
fractal flames (non-linear variations, log-density display, and structural coloring) 
was created in the summer of 1991 while Draves was an intern at the NTT-Data 
Corporation in Tokyo, Japan and was generously allowed to pursue his own 
projects. 

That version was released on the then-nascent world wide web in 1992 under 
the General Public License (GPL), making it the first open-source art. It has 
since been incorporated into and ported to many environments, including the 
Gimp, Photoshop (as Kai’s Power Tools FraxFlame), After Effects, Digital Fu- 
sion, Ultra Fractal 3, screensavers for Macintosh, Windows, and Linux, as well 
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as stand-alone programs (Apophysis, Oxidizer, and Qosmic). 

The combined-channel gamma feature and the vibrancy parameter that con- 
trols it were introduced in 2001. The symmetries were introduced in 2003. Vari- 
ations 7 to 12 were developed by Ronald Hordijk for his screen-saver version of 
the flame algorithm, then ported to the Ultra Fractal version by Erik Reckase, 
and adopted into the original version (with some modifications) by Draves in 
2003. 

In 2004, Mark Townsend released Apophysis, a translation of Draves’ C 
code into Delphi Pascal, with the addition of a GUI for interactive design. Erik 
Reckase became more involved with the flame algorithm after the release of 
Apophysis, and became an official developer and maintainer in 2005. Reckase’s 
contributions have focused on improved image quality, code optimization, and 
keeping up with the wealth of variations being developed in the Apophysis 
community. 

In 2006 Peter Sdobnov added final transforms to Apophysis and they were 
soon after adopted by our implementation to maintain compatibility. 

Thanks to Hector Yee for suggesting tone mapping as the general solution to 
the dynamic range problem, and David Hart for suggesting density estimation 
as an improved filter technique. 

The fractal flame algorithm is also the seed that spawned the Electric Sheep 
distributed screen-saver [2], a follow-on art project by Draves. In this system, 
thousands of idle computers from all over the world are harnessed into rendering 
(and evolving) fractal flames. The work of all participating clients is shared 
alike. 
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Appendix: Catalog of Variations 


For each variation we give its formula and its name, and provide a list of pa- 
rameters for parametric variations. Variables used in these formulas are: 


r= уату 
0 = arctan(x/y) 
Ф = arctan(y/x) 


(а b,c,d,e, f) are considered to be the affine transform coefficients for a 
variation, and are used in dependent variations. Q is a random variable that is 
either 0 or т. A isa random variable that is either -1 ог 1. V is a random variable 
uniformally distributed on the interval [0,1]. The ’trunc’ function returns the 
integer part of a floating-point value. 

A visualization of the distorted coordinate grid and a representative flame 
using this variation are also supplied. As in Figure 2, г and у increase towards 
the lower right, to match the coordinate system of the output image. Note that 
these sample flames are selected based on their characteristic shape and not 
for any aesthetic qualities. The representative images attempt to use a single 
variation, but in general, variations can be mixed when creating flames. For 
random-based variations, a scatterplot is substituted for the coordinate grid. 


Linear (Variation 0) 


Vo(z, у) T (ж, у) 
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Sinusoidal (Variation 1) 


Vi (z, y) = (sinz, sin y) 


Spherical (Variation 2) 


(0,0) = J; (2,9) 
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Swirl (Variation 3) 
V3(x,y) = (z sin(r2) — y cos(r?), z сов(т“) + ysin(r?)) 


Horseshoe (Variation 4) 


Иж) = = ` (z — i) Gg) 2g) 


МА 
N К 
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Polar (Variation 5) 


Handkerchief (Variation 6) 
Volx, y) = т. (sin(0 + r), cos(0 — r)) 


19 


Heart (Variation 7) 
Vz(z,y) = r · (sin(0r), — cos(0r)) 


Disc (Variation 8) 


У: (2,9) = 4 : (sin(zr), сов(лт)) 
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Spiral (Variation 9) 


1 
Vo(z, у) = > cos 0 + sin r, sin 9 — cos r) 


ИИ 


Ж) 


Hyperbolic (Variation 10) 


sin 0 


Vio(z, у) = ( " 


, T COS ) 


КИМ на на 
КУКА 
АД 
КУКА ДО 
МО 


21 


Diamond (Variation 11) 


Vi (x,y) = (sin 0 cos r, cos 0 sin r) 


Ex (Variation 12) 
ро = sin(0 + г), pı = cos(0 — r) 
Vio(x,y) = г> (pg + Pi pò — 91) 
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Julia (Variation 13) 
Via (x, y) = УТ: (сов(8/2 + 0), sin(6/2 + Q) 


(Note: The grid visualization for the julia variation only includes data for 
О = 0.) 


Bent (Variation 14) 


(ж, у) т>0б,у>0 
_ (2х, y) 2а«0у20 
Vule) =) (рујд) 2>0y<0 
( 


2z,y/2 «< 0, < 0 
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Waves (Variation 15) - dependent 


Vis(z,y)— (г ьш (5) ‚у + esin (=)) 


Fisheye (Variation 16) 


Note the reversed order of x and y in the formula. 


2 
Vie(z,y) = —— · (у, z) 
т-ы1 
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Viz(z, y) = (z + csin(tan Зу), у + f sin(tan 32)) 


Popcorn (Variation 17) - dependent 


2 
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Power (Variation 19) 


Vig(a, y) = 7779. (cos 0, sin Ө) 


И 
L 


Cosine (Variation 20) 
Voo(@, y) = (сов(тт) cosh(y), - віп(тт) sinh(y)) 
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Rings (Variation 21) - dependent 
У (2, y) = ((r + с“) mod (2c”) - с“ + r(1 — с?)) · (cos 0, sin 0) 


Fan (Variation 22) - dependent 


+ = пе? 
V _ f r:(cos(0 — t/2),sin(0 — t/2)) (0+ f) mod t > t/2 
3005/8) асов (Оа нін +) (ВЗ T) mod t о 


- 


ESSE 
RA 
NIQO 


ia 
Й 7 
УУЛ ЛУ УУХ, УУ 
ү, ШИН 1 SOUS 


SN 
Aa 


ХОС 
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Blob (Variation 23) - parametric 
pı = blob.high, рэ = blob.low, рз = blob. waves 


ener (» + PEP? (sin(p96) + D) о 


PDJ (Variation 24) - parametric 
pi = pdj.a, p2 = pdj.b, p3 = pdj.c, pa = pdj.d 
Voa(a,y) = (sin(piy) — соя(рәт), sin(paa) — cos(pay)) 
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Fan2 (Variation 25) - parametric 
Fan2 was created as a parametric alternative to Fan. 
pi = 7(fan2.x)”, рэ = fan2.y 


20р2 


D1 


) 


t = 0 + рә — pitrunc( 


· (sin (0 — р1/2) ,сов (0 — р1/2)) t» р1/2 
· (sin (0 + р1/2) ,сов (0 + р1/2)) t € p1/2 


У>5(х, y) = | 
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Rings2 (Variation 26) - parametric 


Rings2 was created as a parametric alternative to Rings. 
p = (rings2.val)? 
t = r — 2ptrunc ЕЗ +7(1—р) 
р 


У>в (1,9) = t: (sin 0, cos 0) 


Еуећзћ (Variation 27) 


Eyefish was created to correct the order of x and у in Fisheye. 


2 
Var (a, у) = d (2 у) 
т+1 
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Bubble (Variation 28) 


Уа (2, у) = ` (ж, у) 


r2 +4 


Cylinder (Variation 29) 


12162 9) = (sina, y) 
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Perspective (Variation 30) - parametric 
р = perspective.angle, рэ = perspective.dist 


p2 
Vao(z,y) = —————À—— - (x, у совр 
30(2, y) зэ (x, у cos pi) 


2 2 H ЫН ius ` ` 
HHHH Toa 
НЫҢ ҒА 
ДАЕТ ААА 
JH ДАДА ДА 
HJT] ААА 
IIFA HAAA 
ААЦ ААА 
ITI] A 
(ЖЕШ ЖЕ LEP NASA 


Noise (Variation 31) 


Vai(z, y) = Фу: (ж сов(2т V5), y sin(21V3)) 
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JuliaN (Variation 32) - parametric 
р = juliaN.power, рэ = juliaN.dist 
рз = trunc(|p;|V) 
t = ($ + 27рз)/рі 


Уз (2,9) = гы. (cos t, sin t) 


JuliaScope (Variation 33) - parametric 
2ı = juliaScope.power, рэ = juliaScope.dist 
рз = trunc([p;|V) 
t = (АФ + 2трз)/р\ 


Узз(х, у) = гы. (cos t, sin t) 
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Blur (Variation 34) 
Vaa(z, у) = Фу. (сов(2т V5), sin(27V5)) 


Gaussian (Variation 35) 


Summing 4 random numbers and subtracting 2 is an attempt at approximating 
a Gaussian distribution. 


4 
Уз(т,у) = | У Wı — 2 | - (cos(27Ws), sin(27Ws)) 
kel 
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RadialBlur (Variation 36) - parametric 


pı = (radialBlur.angle) · (7/2) 
4 


ћ = изв (> Фь- 2), t2 = ó + tı sin pı, із = tı cos pı — 1 
k=l 


1 
Узв (2,9) = s · (r costa + #32, r sin tg + tay) 


Pie (Variation 37) - parametric 


Рі = pie.slices, рэ = pie.rotation, рз = pie.thickness 
tı = trunc(V1p; + 0.5) 


2m 
t2 = р + m + Фәрз) 
1 


Vaz(z, y) = Уз (сов ta, sin 13) 
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Ngon (Variation 38) - parametric 


ngon.circle 


pı = ngon.power, рэ = 27/ngon.sides, рз = ngon.corners, рд 


ta = $ — p2 | %/p2] 


NN 
SS 
CN я 
= = 
^ VI 
со с 
> ~ 
N 
е, 
| 
со со 
> + 
—— 
|| 
bl 
- 


TPL 


Узв (2, y) = К. (x, y) 
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Curl (Variation 39) - parametric 


Рі = curl.cl, рэ = curl.c2 


tı = 1+ piz + р2(2? — 02), t2 = рту + Эрэту 


Узо (2, y) E 


1 
8-0 


: (at, + убо, ytı — ata) 


Rectangles (Variation 40) - parametric 


Рі = rectangles.x, рэ = rectangles.y 


Vao(a, у) = (212/р1 | 


+ 1)pı — z, (Ду/ра] + 1)р2 — v) 
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Arch (Variation 41) 


Tangent (Variation 42) 
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Square (Variation 43) 
Улз(т, у) = (Wi — 0.5, Фэ — 0.5) 


Rays (Variation 44) 


бап(Ф 
Vaa(z,y) = J. (cos x, sin y) 
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Blade (Variation 45) 


Vas(z, y) = т. (сов(Ётодз) + sin(Vrv45), cos(Wrv45) — sin(Wrvas)) 


Secant (Variation 46) 


1 
Ув (т, у) = e ——) 
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Twintrian (Variation 47) 


t = 10810 (sin? (Фтолт)) + cos(Vrva;) 
Vaz(z,y) = x (t,t — wsin(Wrv47)) 


Cross (Variation 48) 


Var(z, у) = VA/(z? — у2)2 - (a, y) 
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